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ABSTRACT 



Aims. We investigate here the X-ray luminosity function (XLF) of absorbed (N/j between 4x 10^' and lO^'' cm ^) and unabsorbed (Nh < 4x 10^' 
^ cm"-) AGN, the fraction of absorbed AGN as a function of Lx (and z), the intrinsic Nh distribution of the AGN population, and the XLF of 
Compton Thick (Nh > 10-"* cm"-) AGN. 

Methods. To carry out this investigation we have used the XMM-Newton Hard Bright Serendipitous Sample (HBSS) a complete sample of 
bright X-ray sources (/,- > 7 x 10"'* erg cm"^ s"') at high galactic latitude (1^1 > 20°) selected in the 4.5-7.5 keV energy band. The HBSS 
. sample is now almost completely identified (97% spectroscopic identifications) and it can be safely used for a statistical investigation. The 
^ HBSS contains 62 AGN out of which 40 are unabsorbed (or marginally absorbed; Nu < 4 x 10^' cm ^) and 22 are absorbed {Nh between 
O I 4 X 10^' and ~ lO^-* cm-^). 

Results. Absorbed and unabsorbed AGN are characterized by two different XLF with the absorbed AGN population being described by a 
steeper XLF, if compared with the unabsorbed ones, at all luminosities. The intrinsic fraction F of absorbed AGN (i.e. the fraction of sources 
with Nh between 4x10^' and 1 O^'' cm"- divided the sources with N/j below 1 0^* cm"- , corrected for the bias due to the photoelectric absorption) 
J^JTj ' with L2-iQkev > 3 X lO''^ erg s"' is 0.57 ± 0.1 1; we find that F decreases with the intrinsic luminosity, and probably, increases with the redshift. 
, Our data are consistent with a flat Log Nh distribution for Nh between 10^" and lO^"* cm"^. Finally, by comparing the results obtained here with 
■ those obtained using an optically selected sample of AGN we derive, in an indirect way, the XLF of Compton Thick AGN; the latter is well 
described by a XLF similar, in shape, to that of absorbed AGN but having a normalization about a factor 2 above. The density ratio between 
Compton Thick AGN (Nh ' 10^-* cm"^) and Compton Thin AGN (Nh ' 10^** cm"^) decreases from 1.08 ± 0.44 at ~ lO^*^ erg s"' to 0.57 ± 0.22 
at ~ lO'*'' erg s"' to 0.23 + 0.15 at ~ 10"*^ erg s"' . 

Conclusions. The results presented here on the anti-correlation between F and -Li are fully consistent with the hypothesis of a reduction of the 
covering factor of the gas as a function of the luminosity and are clearly inconsistent with the simplest unified scheme of AGN. These results 
strongly support the recently proposed radiation-limited dusty torus model although alternative physical models are also consistent with the 
observations. 

Key words, galaxies: active - galaxies: nuclei - galaxies: evolution - X-ray: diffuse background - X-ray: Surveys - X-ray: active galaxies 
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1. Introduction 

For many years the study of the cosmological and statistical 
properties (e.g. luminosity function, cosmological evolution, 
mean spectral properties) of Active Galactic Nuclei ( AGN) was 
principally limited to the optical (see e.g. iBovle et al. 1988; 
Croorn et al. 20041) or soft X-rays (see e.g . Maccacaro et al 
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20051) bands and essentially dealt with complete and unbiased 
samples of optically unabsorbed (type 1) AGN. Indeed the se- 
lection and identification of sources hosting obscured accreting 
nuclei is a difficult task: in the optical domain the active nuclei 
appear very dim and their luminosity could be comparable to 
that of their host galaxies, while in the soft X-ray band (up to 
few keV) their selection is difficult since even hydrogen column 
densities (Nh) of the order of 10^'- 10^^ cm"^ strongly reduce 
the flux from the nuclear source. 
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However, despite their elusiveness, obscured AGN are fun- 
damental for our understanding of the Super-Massive Black 
Holes (SMBHs) history as the large majority of the energy den- 
sity generated by accretion of matter in the Univ erse seem to 
take place in obscured AGN dFabian et al. 1998 ). as testified 
by the integrated energy density contained in the cos mic X- 
ray b ackground (XRB; Setti and Woltjer 1989; Madau et al. 
[19941: IComastriet all 1 19951: iGilli et al.ll2007h . Even more im- 
portant, the recent discovery of quiescent SMBH in the nu- 
clei of no n-active nearby gal axies with prom inent bulges 
dKormendv & Richstond Il995i iMagorrian et al.1 1 19981) . along 
with the presence of scaling relations between the central BH 
mass and galaxy p roperties (e.g. bulge lumi nosity/mass and ve- 
locity dispersion, iFerrarese & MerrittI 1200 0) strongly suggest 
that AGN are leading actors in the formation and evolution of 
gala xies and, in gene ral, of cosmic structures in the Universe 
(see Begeknaiil 2004 and references therein). Missing a large 
fraction of the AGN population could thus bias our understand- 
ing of the evolution of cosmic structure in the Universe. 

Hard X-rays (photon energies between a few keV and ~ 10 
keV) can directly probe AGN activity, since they are almost 
uncontaminated by star formation processes at the X-ray lu- 
minosities of interest (Lx > 10"*^ erg s"' ), and they are sen- 
sitive to absorbed AGN up to an intrinsic absorbing column 
density of A^h~10^"''^"^'* cm , that is to say, they detect all 
but the most absorbed sources. Thus, hard X-ray surveys in 
the 2-10 keV range, as pioneered by AS CA and BeppoSAX 



(Cagnoni e t al. 1998; Ueda et al. 1998; De llaCeca et al.|[l999t 



Fiore et al ] |T999l; lUeda et al. 112003b . and now routinely per- 



formed with Chandra and XMM-Newton^ provide among the 
most complete and unbiased samples of AGN presently pos- 
sible (although they still miss Compton Thick sources i.e. the 
sources with absorbing column densities ! 10^"* cm ). In spite 
of this, a significant fraction of the AGN found in medium and 
deep fields is too faint to provide good X-ray spectral informa- 
tion. Furthermore, the extremely faint magnitudes of a large 
number of optical counterparts of these X-ray faint sources 
make the spectroscopic identifications very difficult, or even 
impossible, with the present day ground-based optical tele- 
scopes. 

Here we use the bright ( fx Z 10"'^ erg cm"^ s"' ), hard (4.5- 
7.5 keV selection band), almost completely identified (97% 
spectroscopic ID) XMM-Newton Hard Bright Sample, to dis- 
cuss a few issues which are currently subject of intense research 
activity, namely: a) the statistical properties (e.g. X-ray lumi- 
nosity functions) of absorbed and unabsorbed AGN; b) the in- 
trinsic ratio of absorbed and unabsorbed AGN as a function of 
the X-ray luminosity; c) the intrinsic absorption column den- 
sity distribution of the AGN population for Hydrogen column 



density up to ~ 10 cm and d) the X-ray luminosity function 
of Compton Thick AGN. These are fundamental issues both to 
study and to follow the accretion history in the Universe and 
to test the unification models of AGN. We stress that the high 
identification rate of the HBSS sample is fundamental to in- 



' A compilation of the main extragalactic surveys executed 
with Chandra and XMM-Newton is reported at the web address 



vestigate these issues since interesting and important classes of 
X-ray emitting sources (e.g. the type 2 QSOs) could be more 
difficult to identify and therefore could be under-represented 
even in samples with an identification rate of the order of 90%. 
Furthermore, good optical and X-ray spectral data are available 
for almost all the sources in the HBSS sample. In particular the 
results of a complete X-ray spectral analysis will be efficiently 
used here to separate absorbed AGN from the unabsorbed ones 
and to compute k-coiTections. 

This paper is organized as follows: in ^we give a short 
account of the XMM-Newton Bright Serendipitous Survey and 
a description of the AGN sample selected in the 4.5-7.5 keV 
range (the HBSS AGN sample hereafter). In ^we study the 
cosmological evolution properties of the HBSS AGN sample 
and derive the X-ray luminosity function (XLF) of the to- 
tal AGN population as well as the XLF of absorbed and un- 
absorbed AGN separately. The intrinsic fraction of absorbed 
AGN as a function of the X-ray luminosities is discussed in ^ 
while in ^we discuss the intrinsic N// distribution of the AGN 
population. In ^we compare our results with the prediction of 
AGN unification models. In ^we compare our results on the 
intrinsic fraction of absorbed AGN as a function of luminos- 
ity with those obtained using an optically selected sample of 
AGN; this comparison will be used to estimate the luminosity 
function of Compton Thick AGN and to derive the ratio be- 
tween Compton Thick and Compton Thin AGN as a function 
of Lx- Finally, summary and conclusions are reported in ^ In 
Appendix A we discuss the method used to take into account 
the photoelectric absorption in the computation of the luminos- 
ity functions. Throughout this paper we consider the cosmolog- 
ical model with (i/o,QM,f^^i)=(65,0.3,0.7); results from other 
papers have been rescaled to this cosmological framework. 

2. The XMM-Newton Bright Serendipitous Survey 
and the IHBSS AGN Sample 

The XMM-Newton Hard Bright Serendipitous Sample used 
here is part of a bigger survey project known as XMM-Newton 
Bright Serendipitous Survey |3(XBS hereafter). This latter con- 
sists of two flux-limited serendipitous samples of X-ray se- 
lected sources at high galactic latitude (\b\ > 20°): the XMM 
BSS sample (389 sources) and the XMM HBSS sample (67 
sources, with 56 sources in common with the BSS sample) hav- 
ing an EPIC MOS2 count rate limit, coiTected for vignetting, of 
10-2 cts/s (2 X 10-^ cts/s) in the 0.5-4.5 keV (4.5-7.5 keV) en- 
ergy band; the flux limit is ~ 7 x 10 '"^ erg cm^^ s"' in boffi 
energy selection bands. At the time of this writing the spectro- 
scopic identification rate is 87% for the BSS sample and 97% 
for the HBSS sample. The details on the XMM-Newton fields 
selection strategy and the source selection criteria of the XMM 
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- The XMM-Newton Bright Serendipitous Survey is one of the 
research programs conducted by the XMM-Newton Survey Science 
Center (SSC, see |http://xmms sc-www.star.le.ac.uk|) a consortium of 
10 international institutions, appointed by ESA to help the SOC in 
developing the software analysis system, to pipeline process all the 
XMM-Newton data, and to exploit the XMM serendipitous detections. 
The Osservatorio Astronomico di Brera is one of the Consortium 
Institutes. 
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BSS a nd HBSS samples are discussed in iDella Ceca et al. 
( 120041) while a description of the optical data and analysis, of 
the optical classification scheme and of the optical properties 
of the extragalac tic sources identified so far is presented in 
Caccianiga et alj (l2008k . The optical ar id X-ray properties of 



the ga lactic population are discussed in lLopez-Santiago et al. 
(l2007h . 



Since the HBSS sample is now almost completely iden- 
tified it will be used in this paper for a statistical investi- 
gation. We note that among the ongoing surveys performed 
with Chandra and XMM-Newton, the HBSS is currently cov- 
ering one of the largest area (~ 25 deg^) and, unlike deep 
pencil beam surveys, is unbiased by problems connected to 
the cosmic variance. The current classification breakdown 
of the HBSS sample is as follows: 62 AGN, 1 cluster of 
galaxies (XBSJ 141830.5+251052) and 2 X-ray emitting stai's 
(XBSJ 014100.6-675328 and XBSJ 123600.7-395217). Two 
X-ray sources (XBSJ 08041 1.3-H650906 and XBSJ 110050.6- 
344331) are still unidentified at the time of this writing. The 
redshift of the 62 AGN are reported in Table 1 along with their 
basic X-ray spectral properties that will be used in this paper 
(e.g. photon index, intrinsic absorbing column density Nh, in- 
trinsic 2-10 keV luminosity). A full account of the X-ray spec- 
tral properties for all the sources in the XMM-Newton Bright 
Survey will be reported in a forthcoming paper. 



2.1. Absorption properties 
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Fig. 1 Intrinsic 2-10 keV luminosity versus intrinsic absorption 
column densities Nh (both quantities derived from a complete 
X-ray spectral analysis) for the AGN belonging to the HBSS 
sample (filled circles: objects for which we have measured the 
Nh', downward arrows: objects for which we have only an up- 
per limit on A^h)- 



The intrinsic 2-10 keV luminosities versus intrinsic absorp- 
tion column densities for the HBSS AGN sample are shown 
in Figure [T] We have directly measured the intrinsic absorb- 
ing column density for 34 AGN while for the remaining 28 
AGN we have only an upper limit. However the measured up- 
per limits are such that we can easily divide the sample into 
absorbed and unabsorbed AGN using a value of Nh = 4 x 10^^ 
cm ~^ . Assuming a Galac tic Ay/Nn ratio of 5.27 x 10"^^ mag 



cm dBohlin et al.lll978h . Nh - 4 x 10^' cm corresponds 



to Ay ~ 2; as shown and discussed in lCaccianiga et al.l (I2008h 
(e.g. see their Figure 6) an Ay value of ~ 2 mag seems to 
be the best dividing value between optically type 1 (Seyl, 
Seyl.5, QSO, NLSyl) and optically type 2 (Seyl. 8, Seyl.9, 
Sey2, QS02) AGN in the XBS survey. Therefore, the choice 
of Nh = 4 X 10^' cm as dividing line between X-ray ab- 
sorbed and unabsorbed AGN has a direct physical link to the 
optical classification. Three objects (XBS J041 108.1-71 1341, 
XBS J205635. 2-0447 17 and XBS J220601. 5-015346) have an 
upper limit between A^h = 4 x 10^' cm and A^// ~ 7 x 10^' 
cm (see Figure [U, so the current X-ray data does not allow 
a firm classification as absorbed or unabsorbed AGN. On the 
other hand these three objects ar e classified as broad lin e (Type 
1) AGN in the optical domain ( Caccianiga et al.ll2008l) : based 
on their optical properties these three objects will be consid- 
ered here as unabsorbed AGN. With these choices the sample 
of unabsorbed (or marginally absorbed) AGN is composed of 
40 objects while the sample of absorbed AGN is composed of 
22 objects. As shown in Section 4 the specific choice of the 



dividing line between absorbed and unabsorbed AGN does not 
affect the main results reported in this paper 

For the sake of clarity we will call in this paper as unab- 
sorbed AGN those AGN with an intrinsic absorbing column 
density A^// below 4x10^' cm"^, as absorbed AGN those objects 
with Nh between 4 x 10^' and ^ 10^"* cm"^ and as Compton 
Thick AGN those sources with Nh above 10^"^ cm"^. 

At the flux limit of the HBSS survey the measured surface 
densities of absorbed (unabsorbed) AGN is 0.87+JJ-23 (1 .59+J^ g) 
deg"^ and the observed fraction of absorbed AGN is 37+7%. If 
we split the sample according to an intrinsic luminosity of 10"^"* 
erg s"' we have that the fraction of absorbed AGN is 29^g^% 
in the high luminosity regime and 42^}^% in the low luminos- 
ity regime, thus confirming first results based on a complete 
sub-sample of 28 HBS S objects fully identified reported in 
Caccianiga et al. I (120041) . However this crude comparison be- 
tween low and high luminosity AGN is biased since it does not 
take into account selection effects due to the absorption, i.e. we 
can observe unobscured sources within a larger volume com- 
pared to the obscured ones since the intrinsic luminosity of the 
latter sources is depressed by photoelectric absorption. We will 
derive and discuss the unbiased intrinsic fraction of absorbed 
AGN as a function of the intrinsic X-ray luminosity in section 
4. 

It is worth noting that the selection strategy of the HBSS 
survey is extremely efficient in finding type 2 QSO, i.e. ab- 
sorbed AGN with an intrinsic luminosity above lO'*'* erg s"' 
. We have 9 of these sources in the HBSS sample (7 with 



4 



R. Delia Ceca et al.: The cosmological properties of AGN in the HBSS 



Nh > 10- 
ulation and 41+[^% 



cm"^); they represent 15^5% of the total AGN pop- 
of the absorbed ones. At the flux limit of 
the HBSS survey the measured surface density of type 2 QSOs 
is 0.36^Q [2 deg"-; this latter surface density of Type 2 QSO 

co mbined with the su rface density of Type 2 QSOs reported 

-2 



in lPerolaetalJ (12004 (~ 48 deg"^ at a 2-10 keV flux Umit of 
) implies a slope in the integral LogN- 

10"'^ and ~ 10"'"* 



10 erg cm ^ s 



LogS of Type 2 QSOs of ~ -2.1 between 



erg cm ^ s ' 



Finally we have not found Compton Thick AGN, which 
are a fundamental population to produce the shape of the XRB 



Gilli etalJ 120071) . 



around its emissivity peak (~ 30 keV, e.g 
As discussed in Appendix A even the 4.5-7.5 keV HBSS sur- 
vey is quite ineflicient in selecting this kind of AGN. Given 
their large flux depression due to the absorption. X-ray surveys 
above 10 keV are nee ded to select this kind of sources (see 
lOella Ceca et al. 2007 for a recent review of this topic). We 
discuss Compton Thick AGN in Section 7 where we derive, in 
an indirect way, their X-ray luminosity function. 



2.2. Redshift and Luminosity distribution 

The redshift and intrinsic 2-10 keV luminosity distribution for 
the AGN belonging to the HBSS sample are reported in Figure 
|2] and Figure [3] respectively. The shaded histogram indicates 
the distribution for absorbed AGN while the empty histogram 
indicates the distribution for the unabsorbed ones. Both classes 
of AGN show a rather flat z distribution with the unabsorbed 
AGN population sampled up to z ~ 1.5 while the absorbed 
AGN population is sampled up to z ~ 0.8. Intrinsic 2-10 keV 
luminosities are distributed over 4 order of magnitude for the 
unabsorbed AGN population (from 10*^ to lO'*^ erg s ' with a 
median logarithmic luminosity of ~ 44.2) and over about 2.5 
order of magnitudes for the absorbed AGN population (from 
~ 7 X 10"^^ to 10"^^ erg s ' with a median logarithmic luminosity 
of 43.7). 



3. Cosmological properties 

3.1. Evolution 




Fig. 2 Redshift distribution of the AGN belonging to the HBSS 
sample (shaded histogram: absorbed AGN; empty histogram: 
unabsorbed AGN) 
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Fig. 3 Intrinsic 2-10 keV luminosity distribution of the AGN 
belonging to the HBSS sample (shaded histogram: absorbed 
AGN; empty histogram: unabsorbed AGN) 



For the sample of unabsorbed AGN (40 obje c ts) we find that, 
accor ding to the Ve/Va test (ISchmidtl[l968l: lAvni & Bahcalll 
the hypothesis of a uniform distribution of the objects 
in flie Universe is rejected at a confidence level of ~ 98% description of the evolutionary prope rties of AGN, bofli in the 
(< Ve/Va >- 0.616 + 0.046). Assuming a pure luminosity 
evolution model with an evolutionary form oc (1 -1- z)*- we find 
a best fit parameter of C^2.7 with an associated 68% confi- 
dence interval of 1.9-3.0. This value of cosmological evolution 
is consistent, within the errors, both with the results obtained 
in the soft (E; 3 keV) energy band using the Einstein Extended 



X-ray energy range (Hasinger e t al.l 120051; lUeda et al. 1 12003 
La Franca et al ]E005 ; Silver man et al .I I2OO7I) and in the optical 
domain dBongiorno et al.ii2007n . 



Medium Sensitivity Survey ( C = 2.56 ± 0.17; lMaccacaro et al 



1991. Delia Ceca et al 



To test this evolutionary behavior, we assume here an 
LDDE model wi th the parametrization as introduced by 
Ueda et al. I ( l2003b . 



1992hand t he EMSS-nRosat AGN sam- 



ples (C - 2.6 ± 0.1; iPage etalJll997.) a nd with the results 
in the 2-10 keV e nergy range reported in lUeda et al. I (120031) 



d<^(L„z) _ d^{L„ 0) 
dLogL-f dLogLx 



X e(z, 



(C = 2.70+°;'^) and lLa Franca etalJ (120051) (C = 3.22+° i^). 



However it is now well established that a Luminosity 
Dependent Density Evolution (LDDE) model provides a better 



"^'^ ^ e(Ze, L.)[(l + z)/(l + Z,)Y\ (Z > Zc) 
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Zc(ix) = 



Z*, (Lx > La) 

zliLxlLaT (Lx < La). 



where Zc corresponds to the redshift where the direc- 
tion of the evolution changes sign. It is worth noting that 
Z(, is a function of the intrinsic luminosity of the object; if 
we assume the best fit parameters of p2=-1.15; z* = 2.49; 
a-0.20; Log L „-45.80 (adapted to Ho=65) as reported in 
La Franca et al] 12005) . then Zc is ~ 0.7, L 1 , L7 for AGN with 
Lx ~ lO'*^, 10**, 10"*^ erg s"' , respectively. Given the coverage 
in the luminosity-redshift plane of the HBSS unabsorbed AGN 
sample, for each luminosity the objects are below Zc, implying 
that we are unable to derive p2 , z*, a and LogLg. For this rea- 
sons we have fixed them from La Franca et al. (2005) and we 
have used the VelVa test to constrain pi. 

We obtain a best fit pi =6.5 with an associated 68% con- 
fidence interval of 3.5 - 10.0. The distribution of the derived 
Ve/Vfl values is consistent with being uniformly distributed 
between and 1 according to a KS test (KS probability ~ 
95%). We have also checked that, given the coverage of the 
luminosity-redshift plane of the HBSS AGN sample, the best fit 
pi is virtually insensitive to the other parameters of the model; 
i.e., pi does not change by varying all the ot her parameters 
within their Icr range as derived from lLa Franca et al. (2005). 

The derived best fit va lue for pi is consisten t, within the 
errors, with that reported in lLa Franca et al.l (120051) (pi =4.62 + 
0.26) and is in very good agree ment with those r e cently ob- 
tained, in the optical domain, by iBongiorno et al.l (l2007b us- 
ing a sample of 130 broad line AGN with redshift up to 
z=5 from the VIMO S-VLT Deep Survey (pi =6.54) and from 
Hopkins et al.l (l2007h using a large data set of AGN selected in 
the Mid-IR, optical, soft X-ray and hard X-ray (pi =5.95+0.23). 

Because of their number statistics (22 objects in total) and 
their distribution in the Lx-z plane the cosmological evolution 
is unconstrained for the absorbed AGN sample (note that the 
absorbed AGN are sampled only up to z~ 0.8). Therefore in 
the following, and in line with the Unification Scheme of AGN, 
we will make the assumption that this class of sources evolve 
with cosmic time (and within the reshift range sampled at the 
HBSS flux limit) in a similar way as the unabsorbed ones. 

3.2. X-ray Luminosity Functions of Absorbed and 
Unabsorbed AGN 

In this section we will derive the intrinsic (i.e. corrected for the 
selection bias due to the intrinsic absorption) local X-ray lumi- 
nosity function of absorbed and unabsorbed AGN. For compar- 
ison with other survey projects the XLF will b e derived in the 



2-lOk eV energy range using the 1/Va method (lAvni & Bahcall 
1980l) as follows. 



For each AGN we have computed the maximum redshift 
(Zma.i) at which the source (characterized by the observed count 
rate in the 4.5-7.5 keV energy band, by the X-ray spectral shape 
and by the redshift) can still be detected above the MOS2 count 
rate limit of the HBSS sample (= 2 x 10"^ cts/s in the 4.5-7.5 
keV band). The X-ray spectral shape of each AGN, as derived 



from the X-ray spectral analysis (reported in Table 1), has been 
used to take properly into account k-corrections. 

This Zmax combined with the solid angle covered by the 
HBSS survey (25.17 deg 2) and the best fit LDDE model de- 
scribed in the previous section allow us to co mpute, for each 
object /, the density-weighted volume (see lAvni & Bahcall 
1980) 



V:. = 25.17 X 



r 

Jo 



e(z, L,) X dV/dz x dz 



The XLFs computed using these volumes are the so called de- 
evolved XLF at z=0. 

In order to obtain the differential 2-10 keV X-ray luminos- 
ity function we bin the individual contributions in bins of equal 
logarithmic width ALogLx according to the intrinsic 2-10 keV 
luminosity as derived from the X-ray spectral analysis of each 
source. For each luminosity bin we have 



dLogLj 



1 



V'ALogL, 



where n is the number of objects in t hat bin. Poisson ian errors 
bars have been computed using the iGehrelsl (119861) prescrip- 
tions. 

In summary, in the XLF computation we have used: 

a) the observed count rate in the selection 4.5-7.5 keV band, 
the count rate limit of the HBSS survey and the best fit LDDE 
model in estimating the available co-moving density weighted 
volume V'a. of each source. The spectral shape of each indi- 
vidual source, as derived from the X-ray spectral analysis, has 
been used to properly take into account k-corrections; 

b) the intrinsic (2-10 keV) luminosities (as derived from the 
X-ray spectral analysis of each source) to divide the sources 
into luminosity bins. 

In Appendix A we discuss a few simulations we have car- 
ried out to prove the validity of this approach. 

In Figure |4] we show (filled dots) a binned representation 
of the de-evolved (2-10 keV) XLF obtained using the total 
AGN HBSS sample (62 objects). The derived XLF (at z=0) 
is compared with the 2-10 keV XLF obtained using local AGN 
samples, the latter being independent from cosmological evo- 
lution. In parti cular in Figure[4 l we co mpare our XLF with that 
obtained from Shinozaki et alj (l2006h using the 2-10 keV se- 



lected AGN in the HEAOl survey ( 49 objects in total, dashed 
line) and by lSazonov & RevnivtsevI (12004) using the 3-20 keV 
selected AGN sample (95 AGN in total; solid line) in the Rossi 
X-ray Timing Explorer (RXTE) slew survey. This last XLF 
has been converted from the 3-20 keV to the 2-10 keV en- 
ergy range using a power-law model with photon index equal 



to 



1.9 ( ^-''°'''' = 0.8). As can be seen in Figure H the best 



fit AGN XLF (at z=0) derived here is in excellent agreement 
with that derived using the RXTE slew data while its normal- 
ization is a factor between 2.5 t o 5 times below that d erived 
using the HEAOl AGN sample dShinozaki et al .1120061) . This 
discrepancy in the norm alization of the local (2- 1 keV ) XLFs 
was already discussed b y lSazonov & Revnivtsev (12004 ') and by 
Shinozaki et al. I (l2006l) and it is currently an open issue. It is 
worth noting that a similar shift between our computation of the 



6 R. Delia Ceca et al.: The cosmological properties of AGN in the HBSS 
Table 2 Best Fit Parameters of the 2-10 keV XLF 



Sample 


Objects 


Log A 


Log 


71 


72 


KS-prob 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Unabsorbed AGN 


40 


-6 21+" '^ 

O.Z.l_Q j2 


44.0 (fixed) 


i.ostl^ 




0.44; 0.24; - 


Absorbed AGN 


22 


-6 62+° " 


44.0 (fixed) 






0.10; 0.49; 0.34 



Columns are as follows: (1) Sample; (2) number of objects in the sample; (3) XLF normalization, A, is in units of h^^ Mpc ^; (4) XLF break 
luminosity in units of h^^ erg s~' . Sinc e is not well constrained w e have fixed it to lO^erg s"' close to that found for the total AGN 
population in lSazonov & RevnivtsevI ( |2004|) (L^, = 10"' erg s"' ) and in lShinozaki et alj ( |2006|) {L^, = 10"'*°^ erg s"' ); (5) low-luminosity 
XLF slope; (6) high-luminosity XLF slope; (7) the three values are the probabilities (ID-KS test) that the observed z, L,- and, in the case of 
absorbed AGN, A''^ distributions are consistent with being derived from the the best fit XLF plus cosmological evolution plus intrinsic 
distribution derived here. See section 3.2 for details. 




Fig. 4 De-evolved (z=0) X-ray luminosity function (2-10 keV) 
obtained using the total HBSS AGN sample (filled circles). We 
have also reported for compa rison the local XLF fo r the total 
AGN population as derived bv lShinozaki et al. (2006) (HEAO 1 
AGN sample; dashed line) and by ISazonov & Revnivtsevi 
(|2004 (RXTE AGN sample; solid line). See section 3.2 for 
details. 



XLF and the IShinozaki et al.l (l2006h XLF is also present when 
we consider only the unabsorbed AGN population, suggesting 
that the possible different distribution of the absorption prop- 
erties of the two samples cannot explain the different observed 
normalizations of the XLFs. Despite the different normaliza- 
tions the shapes of the different XLFs seem to be in satisfactory 
agreement. 

More important is the comparison between the XLF of the 
absorbed and unabsorbed AGN population in the HBSS sam- 
ple reported in Figure |5] (unabsorbed population: open circles; 
absorbed population: filled circles). It is clear that the two pop- 
ulations are characterized by two different XLF with the ab- 
sorbed AGN population being described by a steeper XLF, if 
compared with the unabsorbed ones, at all luminosities. The 



Fig. 5 Comparison between the de-evolved (z=0) X-ray lumi- 
nosity function (2-10 keV) for the unabsorbed (Nh < 4 x 10^' 
cm"^; open circles) and absorbed (4 x 10^' < Nh": 10^* cm"^; 
filled circles) AGN in the HBSS sample. The solid and dashed 
lines represent the best fit two power-law function of the unab- 
sorbed and absorbed AGN XLF, respectively (best fit param- 
eters as reported in Table 2). Please note that for comparison 
reasons the XLF of absorbed AGN has been computed with a 
small shift (A Log L^— 0.03) with respect to that of unabsorbed 
AGN. See section 3.2 for details. 

binned representation of the XLFs reported in Figure |5] have 
been fitted by a smoothly connected two power-laws function 
of the form 

dLogLx Li, L* 

taking into account the error bars of each data point and by min- 
imizing using the routines in the QDP0 software package. 
Best fit XLF parameters and 1 cr errors are reported in Table 2. 



See in http://heasarc.gsfc.nasa.gov/docs/softwarel 
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For the unabsorbed AGN population the value of yi and 
72 are consistent, within the errors, with those reported in 
Sazonov & Revnivtsev (12004') and in IShinozaki et alj (l2006h . 
For the absorbed AGN population the yi (72) derived here is 
slightl y steeper (flatter) than that repo rted in Shinozaki et al. 
feOOe") (ri = l-12!o:!9; n = 3.34+g;gv lSazonov & Revnivtsev 
(2004) does not quote the XLF parameters for the absorbed 
AGN population. 

As discussed in section 2 two HBSS sources 
(XBSJ0804 11.3+650906 and XBSJl 10050.6-344331) are 
still unidentified at the time of this writing. Can their inclusion 
in the AGN sample substantially change the results discussed 
above? In particular their inclusion could be important if these 
sources were high luminosity absorbed AGN where we mea- 
sure a deficit of absorbed sources. To answer to this question 
we have examined in detail their X-ray spectral properties as 
well as the optical (photometric) properties of the most likely 
optical counterpart. XBSJl 10050.6-344331 is well described 
(fixing z=0) by apower-law model having F = 1.80 + 0.16 and 
Nh ^ I +0.5 X 10^1 cm-2. Its 2-10 keV X-ray flux (~ 3 x 10"" 
ergcm"^ s"' ), optical magnitude (mg - 18.0), X-ray to 
optical flux ratio (X/O ~ 1) and X-ray spectral properties 
strongly suggest a type 1 (unabsorbed) object. On the contrary 
XBSJ0804 11.3+650906 is most likely an absorbed AGN. Its 
X-ray spectra is described (at z=0) by a power-law model 
having F = 1.7+0.4 and A^^ = 8.4 ± 4.3 x 10^' cm"^, so with 
an intrinsic Nh (at the source z) well above 4 x 10^' cm"^. The 
optical magnitude of the counterpart (m^ - 21.10) combined 
with its 2-10 keV X-ray flux (~ 2.1 x 10"'^ erg cm"^ s"' 
) implies an X/O flux ratio of ~ 17.5, strongly support ing 



the absorbed AGN hypothesis (e.g. ISevergnini et al. I l2006t) . If 



we use the relationship be tween X/O flux ratio and intrinsic 



luminosity reported in iFiore et al.l (l2003h the most likely 
redshift of this object is ~ 0.55 (consistent with a poor quality 
optical spectrum of the optical counterpart) and its intrinsic 
luminosity is L^^ ^ 1.9 x 10"*"^ erg s"' . The addition of this 
object to the HBSS absorbed AGN sample produce an increase 
of ~ 8% of the space density of the absorbed AGN population 
at Ljc ^ lO'*'* erg s"' , a very marginal diff'erence that is well 
within the error bar reported in Figure |5] It is also worth noting 
that there are four objects (see Figure [TJ which have a best 
fit Nh value below 4 x 10^' cm"^ but with error bars crossing 
this Nh dividing line. The two objects (having large error 
bars) with a luminosity around 2 x 10^^ erg s"' are optically 
classified as broad line QSO; so their optical properties 
are fully consistent with the unabsorbed classification. The 
two objects with a luminosity arou nd 1 x 10"^-^ erg s ~' are 



instead optical elusive AGN dCaccianiga et al.l l2007al) and 



their unabsorbed classification was based on their best fit Nh 
value. On the other end their inclusion in the absorbed AGN 
sub-sample will increase the difference between absorbed and 
unabsorbed AGN around - 10"*^ erg s'^ strengthening 
our conclusions of different XLF of the two populations. To 
conclude, the results presented and discussed in this section 
are stable also taking into account all the uncertainty about the 
two unidentified sources and about our present classification 
(absorbed or unabsorbed AGN) break down. Finally, as 
consistency test, we have compared the observed z, and Nh 



distributions with those predicted by the best fit cosmological 
evolution model derived here (XLF + cosmological evolution) 
folded with the Nh distribution discussed in section 5. The 
ID-KS probabilities reported in Table 2 confirms the goodness 
of our modeling of the AGN cosmological properties. 

4. The intrinsic fraction of absorbed AGN 

One of the most important open issues regarding absorbed 
AGN is to understand their relevance amongst the AGN pop- 
ulation and if (and how) the fraction (F) of absorbed AGN 
changes as a function of Lx and z. We define F as the intrin- 
sic ratio between the AGN with Nh in the range 4 x 10^' - 
lO^"* cm-2 and all the AGN with Nh below lO^'* cm^^. This 
is a l ong standing issue comin g back to almost 26 years ago 
when lLawrence & ElvisI (Il982h reported an anti-correlation be- 
tween the intrinsic luminosity of the AGN and the presence of 
absorption. This issue has been recently re-discussed in many 
papers with contradicting results: several studies su ggest that 



F decreases with L^ and in crease with z (e.g. iLa Franca et al 



20051: iTreister & Urrvl2006h . while other studies suggest that F 



is independent of and z (e.g. iDwelly & P age 200 6i) or fliat 



F is dependent fro r n but not from z ("e.g. lUeda et al. 112003 



Akylas et ai]l2006l loini et al.ll2007i) . Having computed the in 



trinsic de-evolved XLFs of absorbed and unabsorbed AGN we 
can now use them to derive F. We stress that using the de- 
rived XLFs we should be free from the selection bias due to 
the absorption since the appropriate corrections have been al- 
ready taken into account in the computation of the XLF itself. 
We also recall that we are not considering here Compton Thick 
AGN, so the derived F ratio is indeed a lower limit to the true 
fraction of absorbed AGN. Compton Thick AGN will be dis- 
cussed in section 7. 

Using the integral representation of the derived XLFs 
we have first computed F for AGN having an intrinsic lu- 
minosity above ~ 3 x lO'*^ erg s"' (i.e. around the lumi- 
nosity of the faintest obscured AGN in the HBSS sample). 
The fraction (F - 0.57 ± 0.11) is in excellent agreement 
with that found using hard {E > 10 keV) selected (local) 
samples of AGN from INTEGRAL/Swift surveys at a flux 



limit of 



10- 



erg cm 



This comparison is re- 



ported in Table 3, where we have listed the value of F in the 



HBSS sample, in the SW IFT/BAT sample (iMarkwardt et al 



2005 '. lAielloeta 



. ■ l2008h and in the INT EGRAL samples 

( Beckmann et al T l2006t I Bassani et al.l 120061) . Please note that 
the local INTEGRAL/Swift samples are almost free from se- 
lection bias related to the absorption at least for Nh up to 10^^ 
cm"^. The good agreement implies that, after selection effect 
due to the absorption have been taken into account, the HBSS 



is sampling at a fx 



10- 



erg cm s , the same popu- 



lation of absorbed AGN as the present surveys performed at 
£■ > 10 keV. It is also worth noting that after selection effect 
have been taken into account the fraction of obscured AGN in- 
creases from an observed value of 0.37 ± 0.07 (see section 2.1) 
to an intrinsic ones of 0.57 + 0.1 1. 

The differential fraction F of absorbed AGN as a function of 
the intrinsic luminosity is reported in Figure |6]left panel (solid 
line); we have also reported the Icr error bar on F (vertical solid 
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Fig. 6 Left panel: the thick solid hne represents the differential fraction F of absorbed AGN as a function of the intrinsic lumi- 
nosity obtained using the best fit de-evolved XLF of absorbed {Nh between 4 x 10^^ and 10^"^ cm"^) and unabsorbed {Nh below 
4 X 10 cm ) AGN in the HBSS survey. The 68% confidence errors on this ratio are also reported at L, = 10^3, 10'*^, lO^^ 
erg s ' . We have also reported the differential fraction of absorbed AGN obtained assuming no cosmological evolution (short 
dashed line), a cosmological evolution of the absorbed/unabsorbed ratio similar to that recently proposed by iTreister & Urrv 
(12006') (long dashed line) and that obtained using Nh - 10^~ cm ^ as dividing value between absorbed a nd unabsorbed AGN 
(dotted line). The open squa res represent t he fra ction of absorbed AGN as derived from lAkvlas et al. I (l2006h . Middle panel: solid 
line as in the left panel. The Akvlas et al. ( 200d) points have been rescaled to z-0 usi ng the mean redshift in e ach luminosity bin 
and an evolution of the absorbed/unabsorbed ratio similar to that rec ently pro posed bv ITreister & Urrv (l2006h . Right panel: solid 
line as in the left panel. The long dashed line is the prediction from lGilli et al. (2007) on the synthesis model of the CXB. The 
two dashed lines enclose the uncertainty on the F ratio based on the CXB modeling. We stress that all the quantities reported in 
this figure refer to AGN with absorbing column density Nh below 10^"^ cm"-. See section 4 for details. 



Table 3 Fraction of absorbed AGN with L.i ; 3 x 10"*^ erg s ' in 
different samples' 



Sample 




F 


NOTE 


HBSS 


0.57± 


0.11 


XMM 


Markwardt et al. (2005) 


0.59± 


0.09 


SWIFT/BAT 


Beckmann et al. (2006) 


0.54± 


0.10 


INTEGRAL 


Bassani et al. (2006) 


0.56+ 


0.09 


INTEGRAL 


Aiello et al. (2008) 


0.57± 0.13 


SWIFT/BAT 



Notes: ' In deriving, from the original papers, the above reported quan- 
tities we have considered, as for the HBSS AGN sample, only the 
Compton Thin AGN (A'^, < lO^'* cm"^) and a dividing Nh line be- 
tween absorbed and unabsorbed AGN fixe d at 4 x 10^' cm~^. W e have 
not reported the INTEGRAL results from lSazonov et al] ( l2007h since 
these authors do not report the absorbing column densities of their 
objects but simply split the objects into absorbed and unabsorbed ac- 
cording to an Nh value of 10^- cm"-; for comparison the fraction 
of absorbed AGN found from their sample using this Nh value is 
0.42 + 0.09. 



line at - 10"*^, 10'*'*, 10'*^ erg s ' ), computed taking into 
account the current uncertainties on the XLFs. 

Given the coverage in the luminosity-redshift plane of the 
HBSS AGN sample the results reported in Figure |6] are quite 
stable, i.e. they do not depend strongly on the value of the cos- 
mological evolution parameter. To prove this in Figure |6] (left 
panel) we have also reported F as a function of Lx if we assume 
no cosmological evolution for absorbed and unabsorbed AGN 



(short dashed line) or if we assume a cosmological evolution 
of the absorbed population faster than th at the unabsorbed ones 
according to the ITreister & Urrv (l2006h model {F oc {\ + z) , 
long dashed line). Finally we have also reported (dotted line) 
the fraction of absorbed AGN if we re-perform all the analy- 
sis but this time splitting the objects as absorbed or unabsorbed 
using Nh - 10^^ cm"^ as dividing value (44 unabsorbed AGN, 
18 absorbed AGN). 

In the same left panel of Figure |6] we have also re- 
ported the fraction of obscured AGN as recently computed by 
Akvlas et al. I (l2006h using a sample of AGN selected in the hard 
2-8 keV band down to a flux limit of 6 x 10"'^ erg cm"^ s"* 
(open squares), about a factor o ne hundred below ou r flux limit. 
The data points computed by Akvlas et al. (12006 ) have been 



corrected for the different volumes sampled from absorbed and 
unabsorbed AGN (due to the absorption) in a simil ar way as 



done here. Although the HBSS sample and the lAkvlas et al 



(l2006l) show a similar trend of decreasing F as a function of 
the intrinsic luminosity, the fractions computed at z-Q using 
the HBSS sample are systematically bel ow (~ a factor 2 fo r 
Lx > lO'*^ erg s"' ) than those report ed in Akvlas et alJ (l2006h . 
However the data points reported in lAkvlas et all (l2006l)" have 
been computed without taking into account possible differences 
in the cosmological evolution properties of the absorbed and 
unabsorbed AGN population. This effect, that we have already 
shown as of the second order f or the HBSS AGN sample, 
could be very important for the lAkvlas et al. I (12006 ) sample 
since their objects are at a significantly higher z compared to 
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the HBSS AGN sample. To test this effect in Figure |6] (mid- 
dle panel) the lAkvlas et alj ( 12006 ) points have been rescaled to 
z-0 using the mean redshift in each luminosity bins (T. Akylas, 
private communication) and a cosmological evolutio n of F 
according to the model proposed by iTreister & Urrvl (|2006|) 
(F oc (I + z)" "*). A much better agreement with our results is 
clearly evident, suggesting a possible evolution of the fraction 
of absorbed AGN with the redshift. 

Finally in Figure |6] (right panel) we compare F with that 
required to produce the cosmic X-ray background according 
to the new synthesis model reported in lGilli et alj (120071) (long 
dashed line: best fit; short dashed lines: one cr error range; both 
computed and predicted F ratios refer only to absorbed AGN 
with Nh below lO^^ cm^^); we discuss this comparison in sec- 
tion 8. 



5. The intrinsic NH distribution in AGN 

A second important and debated issue on the AGN astrophysics 
is related to their intrinsic Nh distribution, i.e. the Nh distribu- 
tion of the AGN family computed taking into account the selec- 
tion effects related to the photoelectric absorption. To compute 
this distribution we have proceeded as follows. We have first 
split the HBSS AGN sample in four bins of absorbing column 
densities up to A^^ = lO^^ cm^^ (Nh = lO^^ - lO^'; lO^' - lO^^; 
10^^ - lO^"'; 10^-* - 10^"^ cm"^). Second, we have computed the 
integral luminosity function using the objects in each bin of 
Nh and the l/Vg method as discussed in Section 3.2. Third, 
using these XLFs we have determined the density of objects 
in each Nh bin having a luminosity above ~ 3 x 10"*^ erg s"' 
(i.e. around the luminosity of the faintest absorbed AGN in the 
HBSS sample). Finally in order to obtain the Nh distribution 
(reported in Figure [T] upper panel, as solid line) the density in 
each Nh bin has been normalized to the total density of AGN 
with Nh between 10^" and 10^"^ cr n~^. One cr Poisson error 
bars have been computed using the iGehrels prescrip- 
tion and the objects in each bin of A^^- It is now interesting 
to compare the rather flat LogNn distribution obtained from 
the HBSS sample with other samples of AGN or with that re- 
quired from the synthesis model of the CXB. In Figure|7](upper 
panel) we compare the Nh distribution derived using the HBSS 
sample with that computed using hard (£ > 10 keV) selected 
(local) sample of AGN from SWIFT/BAT and INTEGRAL sur- 
veys at a flux limit of ~ 10"" erg cm"^ s"' . We stress again 
that these latter samples are virtually free from selection bias 



related to the absorption for Nh below ^ 10 cm . Clearly 
the flat Nh distribution derived using the HBSS AGN sam- 
ple is in agreement with that derived using local samples of 
X-ray selected AGNat a flux limit about a factor 100 higher, 
thus suggesting that we are indeed sampling a similar popu- 
lation of Compton Thin AGN. The derived flat Log Nh dis- 
tribution is also consistent with th e flat Log Nh distribution 



assumed by iLa Franca et alj (1200 5) in the cosmological anal- 
ysis of the HELLAS2XMM AGN sample. Finally in Figure [T] 
(lower panel) we compare the A^^ d istribution derived here with 



that required by' Gilli et aP ilOQlh (up to Nh = 10^"^ cm"^) to 



fraction of absorbed AGN that is clearly at odd with the results 
reported here. 

It is worth noting that the flat Log Nh distribution der ived 
here is also markely different from that derived by Risaliti e t alj 
(Il999h using a sample of optically selected Seyfert 2 galaxies, 
with the optically selected sample requiring an increasing frac- 
tion of absorbed AGN with increasing the Nh- This difference 
is expected and is due to the fact that we are considering in 
Figure 7 both optically type 1 and optically type 2 AGN, thus 
providing the Nh dist ributi on of the total AGN population. A 
good agreement with iRisaliti et a iJ (119991) is indeed obtained 
is we consider only the HBSS optically Type 2 AGN. In other 



words the Log Nh distribution reported in lRisaliti et al. d 19991) 



does not represent the A^^ distribution of the total AGN popu- 
lation but is probably more typical of that of optically selected 
Seyfert 2 galaxies, which, with few exceptions, are absorbed 
sources with Nh usually above 10^^ cm"^. 

6. Comparison with AGN Unification IVIodels 

The results on the fraction of absorbed AGN as a function of 
the intrinsic luminosity discussed in section 4 can be now di- 
rectly compared with the prediction(s) of the unification model 
of AGN. These results are inconsistent with the simplest unified 
scheme, where a similar thick molecular torus (at a distance of 
the order of one parsec from the central black hole) surrounds 
all the AGN, and which predicts that F is independent from Lx, 
a modifications to this simple zero order unification scheme is 
clearly required. In Figure [8] (upper panel) we compare F-L^ 
relationship derived here with that expected in the framework 
of the original receJw^ torus model (Foe L^*^^, dotted line; 



Lawrencdll991h and in the framework of the radiation-limited 
clumpy dust torus (Foe L dashed line; iHonig & Beckert 
.2007) ; both these models have been normalized with our results 
at Lx - 10'^^ erg s ' .It is clearly evident the excellent agree- 
ment between our results and the prediction of the radiation- 
limited dusty torus discussed in Honig & Beckert (.2007) . 

Finally in Figure [8] (lower panel) we compare the F-L^ re- 
lationship from the HBSS sampl e with the AGN absorption 
model recently discussed bv Lamastra et alj ( l2006i) . According 
to this model the absorption seen in Compton Thin AGN is not 
related to the obscuring torus but to molecular gas in a disk 
located far away from the X-ray source (betw een 25 and 450 



pc dep ending on the BH mass). As shown in iLamastra et al 



(2006) the anti-correlation between the ratio of absorbed AGN 
and the luminosity can be reproduced if a) there is a sufficient 
number of molecular clouds within the radius where the BH 
gravitational influence is dominant and b) there is a statistical 
correlation between the BH mass and the AGN luminosity. We 
have reported in Figure [8] (lower panel) as dashed lines the ex- 
pected behavior of the fraction of absorbed AGN for values 
of the central surface densities of molecular clouds (from bot- 
tom to top) of 40, 80, 120, 200 and 400 Mp/pc^ as reported 
in Figure 4 (right panel) of ILamastra et al. l (l2006h Fl based on 
the assumption Lsoi/LEdd - 0.1 and luminosity dependent- 



fit the spectra of the CXB; the current CXB synthesis modeling 
based on absorbed and unabsorbed AGN require an increasing 



* A. Lamastra lias kindly provided us tiie expected beliavior of their 
model for L^. > 3 x 10'*'* erg s"' . 
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Fig. 7 Upper panel: intrinsic Nh distribution between 10 and 
10^"^ cm"^ obtained using the HBSS AGN sample (see sec- 
tion 5 for details), c ompared with Nh dis tribution from the 
SWIFT/BAT sample jMarkwardt et al. l2005t long dashed line' 



and from the INT EGRAL samples ([Beckmann et al.l 1200^ 



short dashed Une; iBassani et al.l l2006t dotted line). Lower 



panel: intrinsic Nh distribution between 10^' and 10^^ cm ^ 
obtained using the HBSS AGN sample, compared wit h that 



requir ed from the synthesis mo del of the CXB (Gilli et al 



(120071) ) ■ In this panel, as done in lGilli et al.ll2007L the number 
of objects in each Nh bin is normalized to the total number of 
AGN with Nh between 10^' and 10^^ cm^^. 



bolometric corrections (see iLamastra et al. ]|2006' for details). 
Our results on the F-Lj^ relationship are fully consistent with 
this model if the surface density of molecular clouds in the cen- 
tral part of the galaxy is betwee n 80 and 2 00 Mq / pc^, a con- 



dition that, as discussed in Lamastra et al.l 62006.) is consistent 
with the observations of local galaxies. 

7. The X-ray luminosity function of Compton Thick 
AGN 

As discussed above (and as done so far in this paper) XMM- 
Newton and Chandra data can be efficiently used to investigate 
the statistical properties of AGN with column densities below 
~ lO^** cm"^. On the contrary the statistical properties (e.g. 
XLF) of absorbed sources having Nh above ~ lO^'* cm"^, the so 



Fig. 8 Upper panel: the thick solid line represents the fraction 
F of absorbed AGN as a function of the intrinsic luminosity 
obtained using the best fit de-evolved XLF of absorbed (Nh 
between 4 x 10^' and 10^"* cm"^) and unabsorbed {Nh below 
4 X 10^' cm 2) AGN in the HBSS sample. The 68% confidence 
errors on this ratio are also reported at - lO"*^, lO''^, 10*^ 
erg s ' . The dotted line represents the absorbed fraction ex- 
pected from the original receding torus model discussed in 
Lawrence! (Il99lh {F oc L "^), while the dashed line repre- 
sents the absorbed fraction expe cted from the radiation-li mited 
dusty torus recently discussed in lHonig & BeckertI (l2007h {F oc 
^-0.25-). fjjese expected fractions have been normalized to the 
observed values at - lO'*-' erg s"' ; Lower panel: the thick 
solid line as in the upper panel. The dashed lines represent the 
expected behavior of the fraction of absorbed AGN for values 
of the galaxy central surface densities (from bottom to top) of 
40, 80, 1 20, 200 and 400 M^/p c^ according to the models dis- 
cussed in lLamastra et al. ( 20061) . 



called Compton Thick AGN, are currently completely uncon- 
strained. Acc ording to the late st versions of synthesis modeling 
of the XRB (iGilli et al.ll2007l) . these sources should represent 
a substantial fraction of the total AGN population but only a 
few dozens of them (mostly local) have been found and stud- 
ied so far in the X-ray domain (see Delia Ceca et al.ll2007 for a 
review). 

In Figure |9] we show the fraction of optically narrow line 
AGN as a function of the X-ray luminosity (dashed line) as de- 
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43 44 

log L, (2-10 keV; erg s"') 



Fig. 9 The dashed line represents the fraction of optic ally 
narrow line AGN {F optical) as derived by SimpsonI (12005 ) us- 
ing a complete, magnitude-limited, sample of active galax- 
ies from the Sloan Digital Sky Survey and assuming the 



mean Li 



-IQkev 



ratio for Seyfert galaxies derived by 
Mulchaev et alj (119941) . The thick solid line represents the frac- 
tion F of absorbed AGN (sources with Nh between 4x10^' and 
10^"* cm"^ divided by the sources with Nh below lO^'* cm"^) as 
a function of the intrinsic luminosity obtained using the best 
fit de-evolved XLF of absorbed and unabsorbed AGN in the 
HBSS sample. The 68% confidence errors on this ratio are also 
reported at = 10'*\ lO"*"*, lO'*^ erg s"' . 



rived by SimpsonI ( 2005f ) using a complete, magnitude-limited, 
sample of active galaxies from the Sloan Digital Sky Survey 
and assuming the mean Z^r p///i/^(2-iofeiO - 0.01 5 ratio for 
Seyfert galaxies obtained by Mulchaev et alj ( 19941) (fully con- 



sistent with the similar value reported in iHeckman et alj|2005 



for the unobscured view of Seyfert galaxies, L[oni] /^(i-iofev) - 
0.017). The thick solid line (and the error bars) are the re- 
sults derived here on the fraction of absorbed (4 x 10^^ < Nh 
1 10^4 cm 2) AGN obtained using the HBSS AGN sample (see 
Section 4); we recall here that in the HBSS we are sampling 
only AGN with Nh below lO^'* cm"^. 

The be havior of the frac tion of optically narrow line AGN 
derived by 'Simpsonj ( 2005 ) is remarkable similar to the cor- 
rela tion betw e en F a nd Lx found here. More important, since 
the SimpsonI (2005) sample should contain both Compton 



Thin (Nh 1 lO^^ cm '^) and Compton T Tiick (Nh 1 10^^ cm 



,-2 
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AGN (see also the results discussed in IHeckman et al. I2005i) . 



soni 



the comparison between the results found by ISimpsonI (120051) 
(dashed line in Figure|9j ¥ optical hereafter) and those found here 
(solid line in Figure |9| allow us to estimate the ratio between 
Compton Thick AGN and AGN with Nh between 4x10^' and 



10^"^ cm as a function of the luminosity. 



,-2 



Assuming that the number of Compton Thick AGN is 
proportional to the number of AGN with Nh in the range 
[4 X 10^'; ~ 10^"^ cmr\ Nnick - C x Naeix-eia, thus from the 
comparision between F and Fopticai, and recalling that the defi- 
nition of F given in Section 4 does not include Compton thick 



AGN in either the numerator or denominator, we can derive 

(Fopticai ~ F) 



c 



Fx (I -Fopticai) 



In order to derive an estimate of the density of Compton 
Thick AGN as a function of the luminosity we have first com- 
puted C for three values of the intrinsic X-ray luminosity 
(L,. = 10"*^ L, = lO'^'* and = lO'*^ erg s"' ); thus we have 
multiplied them for the corresponding XLF of absorbed AGN 
with Nh between 4 x 10^' and ~ 10^"* cm"^ (i.e. second line of 
Table 2). The obtained results are reported in Figure[TO]as open 
squares (the errors have been computed taking into account the 
errors reported in Figure |9ll. 

These estimates of the density of Compton Thick AGN at 
L, = 10'^^ L, = lO'** and L,. = 10"*^ erg s"' can be well 
reproduced assuming an XLF similar, in shape, to that of ab- 
sorbed AGN reported in Table 2 but having a normalization 
Arhick = 2 X A4£2i-£24; this XLF is shown as a thick solid 
line in Figure [10] while the two dashed lines corresponds to 
the XLFs with Arhick = A4E21-E24 (lower dashed line) and 
Arhick = 4 X A4E2\-E24 (upper dashed line). 

In Figure [To] we have also reported: 

a) the density of local Compton Thick AGN (filled circles 
at Lx ~ 3 X IQ*^ and -3x10^"^ erg s"' ) as derived using the 
three Compton Thick AGN (NGC 4945, NGC 3 281 and MKN 



3) in the INTEGRAL all sky survey at \b\ > 5° (ISazonov et al 
2007h To convert the observed 17-60 keV fluxes into 2-10 



keV intrinsic (i.e. unabsorbed) luminosities we have corrected 
the fluxes for intrinsic absorption and assumed F = 1 .9; 

b) the density of Compton Thick AGN with Lx ~ (1 - 
4) X 1 0"^^ erg s ' (filled triangle) as derived by iDaddi et al 



(120071) using Spitzer and Chandra data in two GOODS fields 
with multiwavelenghts information. Their estimated space den- 
sity at z between 1 .4 and 2.5 has been rescaled to z=0 using our 
best fit evolutionary model; 

c) the density of Compton Thick AGN (filled squares) with 
Lx ~ a - 10) X 10"*^ ergs-' (z between 0.7 and 1.0) and 
Lx ~ (1 - 10) X 10^^ erg s"' (z between 1.3 and 2.0) as de- 
rived by lFiore et al. (2008) using the Spitzer and Chandra data 
in the COSMOS field. Their estimated space density has been 
rescaled to z=0 using our best fit evolutionary model. 

As it is clear in Figure [TO] the different estimates of the 
space density of Compton Thick AGN, derived using differ- 
ent seletion criteria (and methods), are in remarkably good 
agreement with the suggested XLF of Compton Thick AGN 
(n = 1.55, 72 = 2.61, Log = 44, Log Arhick = -6.32) wifli a 
possible hint of a XLF flattening below few times lO'*^ erg s"' 
(where, we point out, the XLF of absorbed AGN in the HBSS 
is not well defined since we have no objects in the sample). 

The good agreement seem to suggest that, at the sam- 
pled X-ray luminosities, the space density of optically elu- 
sive Compton Thick AGN (e.g. objects like Arp 299, see 
Delia Ceca et al. 12002") s, at most, similar to that of optically 
selected Compton Thick AGN since the former population 



^ We thank S. Sazonov to have provided us the sky coverage of this 
survey in a tabular form. 
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43 44 45 

Log (2-10 keV; erg s"') 

Fig. 10 The open squares represent the density of Compton 
Thick AGN at Lx = 10'*\ lO^^ and 10"*^ erg s"' derived as 
explained in section 7. The thick solid line represents the pro- 
posed XLF of Compton Thick AGN with A r/iici: = 2xA4£21-£24 
(n = 1.55, 72 = 2.61, Log L* = 44, Log A™^ = -6.32) while 
the two dashed lines corresponds to the XLFs with Aruck = 
A4E21-E24 (lower dashed Une) and Arhwk = 4 x A4E21-E24 (up- 
per dashed line). The filled circles at L(2-iokeV) ~ 3 x lO'*^ and 
~ 3x10^"^ erg s"' represent the density of local Compton 
Thick AGN derived using the Compton Thick AGN in the 
INTEGRAL all sky survey nl\b\ > 5°. The filled ti-iangle rep- 
resents the density of Compton Thick AGN with Lx ~ (1 - 



4) X lO'*^ erg s ' as derived bv lDaddi et all (l2007h . The filled 



squares represent the density of Compton Thick AGN with 
Lx ~ (3- 10)xlO' 



derived by lFiore 



Oferg_s_iand Lx ~ (1 - 10) x lO'*^ erg s"' as 
et alJ(l2008l). The den sities of Compto n Thick 



AGN reported in 'Daddi et^ (120071) and iFiore et all (120081) 



have been rescaled to z-0 using our best fit evolutionary model. 
See section 7 for details. 



should be clearly missi ng in the optical ly selected sample of ac- 
tive galaxies defined in lSimpsonI (l2005i) but should be present 



in the infrared or X-ray selected samples. This fact is also sup- 
ported by the consideration that the total number of Compton 
Thick AGN cannot be increased arbitrarely but we have to 
take into account the limits imposed by the local black mass 
density derived by lMarconi et al.l (120041) from dynamical stud- 
ies of local galaxy bulges (pbh - ^-^Uhl^ X 10^ Mfr, Mpc '-^ 



for Hn =65). Using the formalism described in lLa Franca et al 
(l2005h (see their equations 15, 16 and 17), assuming a radia- 
tive efficiency ~0J_, a bolometric conversion factor equal to 
25 dPozzi et alJl2007l) and integrating the XLFs of AGN (un- 
absorbed, absorbed and Compton Thick) from z=4.5 to z=0 
and from Lx ~ 3 x 10"^^ erg s"' to ~ 10"*^ erg s"' , the limit 
of pBH - 5.6/!g5 X 10^ M© Mpc"^ is violated if the XLF of 



Compton Thick AGN is more the ~ 4 times that of absorbed 
AGN, e.g. the upper envelope reported in Figure 10; we stress 
tha this is probably an upper envelope since we have not con- 
sidered in the computation of the local black mass density the 
objects with luminosity below ~ 3 x 10"*^ erg s"' . 

Finally the comparison of Foptkai and F reported in Figure|9] 
can be used to evaluated the ratio, Q, between Compton Thick 
and Compton Thin {Nh < 10^"* cm"^) AGN. Assuming, as 
above, that the number of Compton Thick AGN is proportional 
to the number of AGN with Nh in the range [4 x 10^' ; ~ lO^'* 
cm"^], thus 



Q 



(f, 



optical 



-F) 



(1 Fopttcai) 



We derive that the density ratio between Compton Thick 
AGN and Compton Thin AGN decreases from Q = 1 .08 + 0.44 



at 

Q 



. io« 

0.23 d 



erg s to 2 = 
0.15 at ~ lO'*^ 



0.57 ± 0.22 at 
erg s"' . 



10^^ erg S-' to 



8. Summary and Conclusions 

We have discussed here the cosmological properties of the 
62 AGN belonging to the XMM-Hard Bright Serendipitous 
Survey, a complete and representative sample of bright ( fx > 
7 X 10"'"* erg cm"^ s"') serendipitous XMM-Newton sources 
selected in the 4.5 - 7.5 keV energy band on a sky area of ~ 25 
deg"^. Since the HBSS sample is almost completely spectro- 
scopically identified (ID rate ~ 97%) it allows us to have an 
unprecedented and unbiased view of the extragalactic 4.5 - 7.5 
keV sky in the bright flux regime. Using an Nh dividing value 
of 4 X 10^' cnT^ (see Section 2.1), the HBSS AGN sample 
is composed of 40 unabsorbed (or marginally absorbed) AGN 
and 22 absorbed AGN. 

The main results reported and discussed in this paper are: 

1. the HBSS survey is extremely efficient in finding type 2 
QSO, i.e. absorbed AGN with an intrinsic luminosity in ex- 
cess to lO"*"* erg s"' . They represent (15 + 6)% of the total 
AGN population and (41 ± 13)% of the absorbed ones. At 
the flux limit of the HBSS survey (~ 7x10"'"* erg cm"^ s"' 
in the 4.5-7.5 keV energy range) the measured surface den- 
sity of type 2 QSOs is 0.36+°;j^ deg-^ (0.28+°;j^, deg-^ if 
we consider only the 7 Type 2 QSO with Nh > 10^^ cm"^). 
At the same flux limit the measured surface density of ab- 
sorbed (unabsorbed) AGN is 0.87+[;-23 (1.59+Jj |) deg-^ and 
the observed fraction of absorbed AGN is (37 + 7)%; 

2. the de-evolved (z=0) 2-10 keV X-ray luminosity function 
of the total AGN population derived using the AGN HBSS 
sample is in very good agreement with the local AGN XLF 
obtained by ISazonov & RevnivtsevI (120041) using the data 
from the Rossi X-ray Timing Explorer survey. A smaller 
normalization (a factor between 2.5 to 5 times depending 
on L^) is observed if compared with the local HEAOl A GN 
XLF as recently computed from lShinozaki et al. I (l2006h : 

3. absorbed and unabsorbed AGN are characterized by two 
different XLF with the absorbed AGN population being de- 
scribed by a steeper XLF, if compared with the unabsorbed 
ones, at all luminosities; 
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4. the intrinsic fraction of absorbed AGN with Lx > 3 x 10*^ 
erg s"' (i.e. around the luminosity of the faintest absorbed 
AGN in the HBSS sample) is 0.57 + 0.11. This is in excel- 
lent agreement with that found using hard (£ > 10 keV) se- 
lected (local) sample of AGN from INTEGRAL/Swift sur- 
veys at a flux limit of 10"" erg cm~^ s"' ; 

5. the fraction F of absorbed AGN is a function of lumi- 
no sity. A comp a rison of our results with those reported 



Akvlas et al. (l2006h (obtained using a source sample 



at fainter fluxes and thus at higher z) seem to suggests 
an evolution of the fraction of absorbed AGN with z. 
This evolution is consistent with that recently proposed by 
llreister & Urry (2006); 

6. we find a flat Log Nh distribution for Nh between 10^" and 
10^"* cm"^, still in very good agreement with local samples 
of hard (£ > 10 keV) selected AGN at a flux Hmit of 10 " 
erg cm"^ s"' . 

7. the shape of the F-Lj relationship and the derived flat Log 
Nh distribution (both in go od agreement with the assump- 
tions and the results o f|La Franca et alj|2005h are not con- 
sistent with the current models for the production of the 
cosmic X-ray background based on the absorbed AGN. 
This disagreement, probably due to the fact that we are 
not sampling the X-ray source population producing the 
bulk of the CXB, may require a revision in the current 
model's parameters, e.g. an increase of the fraction of ab- 
sorbed sources with redshift and/or a Nh distribution which 
depends on the source luminosity, redshift or both; 

8. by comparing the results obtained here with the behav- 
ior of the fraction of optical narrow line AGN derived by 
[Simpson (2005) using a complete, magnitude-limited, sam- 
ple of active galaxies from the Sloan Digital Sky Survey we 
have derived, in an indirect way, the X-ray luminosity func- 
tion of Compton Thick AGN; it is described by a smoothly 
connected two power-law function (see section 3.2) having 
n = 1.55, 72 =2.61, Log = 44, Log Amck = -6.32. The 
density ratio, Q, between Compton Thick AGN (A^^! 10^"* 
cm"^) and Compton Thin AGN {NhI 10^"* cm"^) decreases 
from 2 = 1.08 + 0.44 at ~ 10"*^ erg s"' to g = 0.57 + 0.22 
at ~ lO'*'* erg s"' to g = 0.23 + 0.15 at ~ 10"*^ erg s^^ ; 

9. finally, the decreasing fraction of absorbed AGN as a func- 
tion of the luminosity is fully consistent with the hypoth- 
esis of a reduction of the covering factor of the gas com- 
ponent as a function of the luminosity. A similar reduction 
with the luminosity of the covering factor of the circum- 
nuclear dus t component has been recently pointed out by 



Maiolino et al.l (l2007l) using a sample of AGN observed 



with Spitzer. These results are clearly inconsistent with 
the simplest unified scheme of AGN. Indeed we found 
an excellent agreement between our results and the pre- 
dictions o f the radiation-limited dusty torus discussed in 



Honig & B eckert (2007). An alternative model which also 



explains well our re sults has been recently discussed by 
Lamastra et alj ( 1200 6) and is based on the obscuration due 



to molecular gas in a disk located far away from the X-ray 
source (between 25 and 450 pc depending on the BH mass). 
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Table 1: Basic X-ray information of the XMM-Newton HBSS AGN 
Sample 



Source 



(1) 



Rate 

cts/s 
(2) 



(3) (4) 



Nh 
xl022 

-2 

cm 
(5) 



+mn 

-0.06 

(frozen) 

+0.14 



m 

).Q4 
).04 

m 
i? 

]'.ol 



fl-WkeV 

erg cm"^ s"' 
(6) 



Log L2-\0keV 

erg s~' 
(7) 



XBSJ002618 
XBSJ013240 
XBSJ013944 
XBSJ015957 
XBSJ021640 
XBSJ021808 
XBSJ021817 
XBSJ021822 
XBSJ023713 
XBSJ030206 
XBSJ030614 
XBSJ031015 
XBSJ031146 
XBSJ031859 
XBSJ033845 
XBSJ040658 
XBSJ040758 
XBSJ041108 
XBSJ050536 
XBSJ052108 
XBSJ052128 
XBSJ074202 
XBSJ074312 
XBSJ083737 
XBSJ083737 
XBSJ091828 
XBSJ095218 
XBSJ101850 
XBSJ101922 
XBSJ104026 
XBSJ 104522 
XBSJ104912 
XBSJ112026 
XBSJl 13106 
XBSJ113121 
XBSJ113148 
XBSJ122656 
XBSJ124641 
XBSJ132038 
XBSJl 33942 
XBSJ 134656 
XBSJl 34749 
XBSJ140102 
XBSJ140113 
XBSJ140127 
XBSJ141531 
XBSJ 142741 
XBSJ143835 
XBSJ143911 
XBSJ153452 



5+105019 
.1-133307 
.0-674909 
5+003309 
.7-044404 
,3-045845 

4- 045113 
2-050615 

5- 522734 
8-000121 
1-284019 

5- 765131 

1- 550702 

2- 441627 

7- 352253 

8- 712457 

9- 712833 
1-711341 

6- 290050 

5- 251913 
9-253032 
7+742625 
1+742937 
0+255151 
1+254751 
4+513931 
9-013643 
5+411506 
6+412049 
9+204542 
1-012843 
8+330459 
7+431520 
9+312518 
8+310252 
7+311358 
5+013126 
8+022412 
0+341124 

6- 315004 
7+580315 
9+582111 
0-111224 
4+024016 
7+025605 
5+113156 
8+423335 
1+642928 
2+640526 
3+013104 



2.34 
3.23 
2.05 
3.80 
2.08 
2.67 
3.30 
4.54 
3.23 
3.10 
4.61 
4.39 
5.87 
2.16 
2.37 
3.41 
4.96 
2.20 
2.19 
2.11 
3.08 
3.38 

10.92 
2.98 
7.30 
3.03 

24.50 
2.16 
2.46 
8.70 
3.00 
2.06 
2.32 
2.27 
3.88 
3.43 
3.07 
2.22 
2.89 
3.52 
3.33 
7.39 
7.21 
2.10 
6.66 
2.58 
3.62 
3.57 
2.44 
8.36 



+ 0.55 
±0.71 
±0.46 
+ 0.92 
±0.49 
±0.25 
+ 0.40 
±0.43 
±0.61 
±0.42 
±0.92 
±0.46 
± 1.12 
±0.49 
±0.36 
±0.68 
±0.91 
±0.42 
±0.44 
±0.55 
±0.88 
±0.49 
±0.64 
±0.71 
± 1.12 
±0.54 
±2.98 
±0.49 
±0.47 
± 1.15 
±0.67 
±0.51 
±0.37 
±0.44 
±0.69 
±0.48 
+ 0.54 
±0.49 
±0.40 
±0.51 
±0.56 
±0.66 
±0.59 
±0.52 
±0.62 
±0.34 
±0.53 
±0.54 
±0.39 
± 1.30 



0.473 
0.562 
0.104 
0.310 
0.873 
0.712 
1.080 
0.044 
0.193 
0.641 
0.278 
1.187 
0.162 
0.140 
0.113 
0.181 
0.134 
0.923 
0.577 
1.196 
0.588 
0.599 
0.312 
0.105 
0.080 
0.185 
0.020 
0.577 
0.239 
0.465 
0.782 
0.226 
0.146 
1.482 
0.190 
0.500 
0.733 
0.934 
0.065 
0.114 
0.373 
0.646 
0.037 
0.631 
0.265 
0.257 
0.142 
0.118 
0.113 
1.435 



2.02 
1.90 
1.94 
2.14 
2.24 
2.01 
1.84 
1.66 
1.90 
1.90 
1.56 
1.92 
2.08 
1.72 
1.90 
1.90 
1.90 
1.95 
1.85 
1.72 
1.90 
2.02 
1.99 
1.77 
1.92 
1.47 
1.90 
2.29 
2.13 
1.99 
2.05 
1.68 
1.65 
1.77 
1.62 
1.90 
1.61 
2.22 
1.77 
1.62 
1.90 
2.19 
1.90 
2.05 
1.57 
1.84 
1.90 
1.84 
1.90 
1.75 



-0.04 
+0.13 
-0.12 
+0.43 
-0.37 

(frozen) 

(frozen) 
(frozen) 

+0.48 



-0.19 

(frozen) 

+0.16 

+8:oi 

-M 
-M 

-0.12 
+0.51 
-0.30 

(frozen) 

+0.09 

-0.10 
+0.05 

-M 
+8:8? 

-0.42 
+0.37 

-M 

-0.28 

(frozen) 

+0.23 

-M 

-0.18 

(frozen) 

+0.02 



-0.24 
+0.05 



-0.05 

(frozen) 

+0.42 
-0.18 

(frozen) 

+0.09 

=QM 



< 0.011 
2 540 

< 0.018 

< 0.010 

< 0.020 

< 0.007 

< 0.040 
20 540 +•^•2*° 

< 0.180 

< 0.015 

< 0.048 

< 0.022 

< 0.013 

390 +"■3'*° 

yj.jyyj _o.280 

28.500 IiIqqo 
21.970 !J5| 
28.230 

< 0.610 



0.610 



+0.210 



0.150 



100 
12 710 
070 

< 0.040 
340 

m 

-0Q4Q 
4 740 

29 400 

< 0.025 

< 0.040 

< 0.009 

< 0.038 
0.030 r"050 
6.280 : 
0.065 



+o:S8 

760 ^"-^^^ 
/uu _Q 30Q 

3 200 '^'^■^'^ 

J.Z,L>U _g46g 

< 0.038 

150 

230 ^"-^ 
9 500 

< 0.015 

< 0.004 

< 0.170 
150 

U.IJW _Q Q20 

< 0.022 
4.480 :^°"30 
1.850 



-0,550 

20 000 ™ 
< 0.050 



2.83 
1.76 
1.15 
2.91 
1.11 
2.32 
2.65 
3.00 
2.74 
2.23 
2.91 
3.41 
2.79 
1.63 
1.67 
1.59 
2.56 
0.82 
1.33 
2.46 
1.34 
1.64 
9.67 
3.04 
6.63 
2.62 
11.00 
1.51 
2.87 
6.21 
2.11 
1.83 
1.87 
1.08 
2.70 
2.96 
2.23 
1.36 
2.54 
1.70 
1.47 
5.14 
5.21 
0.34 
7.42 
2.09 
2.11 
2.33 
1.09 
7.56 



44.44 
44.43 
42.56 
44.03 
44.74 
44.79 
45.22 
42.53 
43.52 
44.63 
43.87 
45.47 
43.37 
42.99 
43.21 
43.56 
43.44 
44.59 
44.27 
45.26 
44.45 
44.46 
44.55 
43.00 
43.09 
43.52 
42.54 
44.43 
43.71 
44.76 
44.85 
43.48 
43.20 
45.14 
43.50 
44.54 
44.73 
44.90 
42.48 
42.81 
44.04 
45.07 
42.28 
43.83 
44.23 
43.68 
43.20 
43.03 
42.96 
45.95 
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Table 1: continued. 



Source 








Rate 


z 


r 




Nh 


ii- 


-lOkeV 


Log L2-mev 










xlO-3 








Xl022 


xlO-13 












cts/s 








cm"^ 


erg cm' 




erg s"' 


(1) 








(2) 


(3) 


(4) 




(5) 




(6) 


(7) 


XBSJ160645 

XBSJ161820 


.9+081525 

.7+124116 


7, 

2, 


,36 

,09 


± 1.37 
+ 0.61 


0.618 

0.361 


1.75 

1.90 


+0.71 
-0.71 

(frozen) 


14 270 

4 960 

^.yvu 2.480 




j.o4 

0.88 


44.94 
43.72 


VD CT1 i^C/lOC 


.3 + 142159 


c 

J, 


,26 


1 1 AO 

± l.Oo 


A 1 no 

0.1 /o 


oil 
2.1 1 


+0.07 
-0.03 


< 0.016 




6.34 


43.63 


XBSJ165448 


.5+141311 


6, 


,25 


± 1.79 


0.320 


1.84 


+0.13 
-0.08 


< 0.043 




4.46 


44.22 


XBSJ193248 


.8-723355 


4, 


,47 


+ 0.73 


0.287 


1.48 


+0.25 
-0.22 


730 

u. / JU .0 290 




2.30 


43.80 


XBSJ204043, 


.4-004548 


3, 


,24 


+ 0.80 


0.615 


1.90 


(frozen) 


3 000 

J.\J\J\J _o.900 




1.66 


44.50 


XBSJ205635, 


.7-044717 


2, 


,08 


±0.50 


0.217 


1.91 


+0.51 

M 


< 0.560 




1.66 


43.42 


XBSJ205829, 


.9-423634 


3, 


,91 


±0.56 


0.232 


1.91 


0.090 

0.080!: 

510 




3.17 


43.76 


XBSJ2 13002, 


.3-153414 


2, 


,30 


±0.47 


0.562 


2.10 






1.84 


44.45 


XBSJ213820, 


,2-142536 


2, 


,80 


±0.59 


0.369 


1.61 


m 




2.26 


44.05 


XBSJ214041, 


,4-234720 


3, 


,30 


±0.68 


0.490 


2.19 


m 


< 0.012 




1.69 


44.28 


XBSJ220601, 


,5-015346 


2, 


,26 


±0.55 


0.211 


1.56 


-0.09 


< 0.670 




1.66 


43.36 



Columns are as follows: (1) Source name; (2) Source MOS2 count rate, and Icr error, in the 4.5-7.5 keV energy band (units 
of 10"^ cts/s). Please note that the reported count rates have been corrected for vignetting; (3) Redshift; (4) Source photon index, 
and 90% confidence errors; (5) Absorbing column density at the source redshift, and 90% confidence errors; (6) Observed 2-10 
keV flux corrected for Galactic absorption; (7) Intrinsic (i.e. de-absorbed) 2-10 keV luminosity. 
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Appendix A: Computation of the intrinsic 2-10 keV 
XLF using the IHBSS AGN sample 

In Section 3 we have derived the intrinsic 2-10 keV XLF of 
absorbed and unabsorbed AGN using the HBSS AGN sam- 
ple selected in the 4.5-7.5 keV energy band. To this purpose 
we have used the 1/Va method, a well established procedure 
to derive space densities of interesting astronomical objects 
dAvni & Bahcall 198 0). 

However to apply this method to the absorbed AGN popu- 
lation the following considerations have to be taken into ac- 
count. Given an intrinsic count rate (i.e. the count rate we 
would have observed without intrinsic absorption), the ob- 
served count rate is a function of the absorbing column density 
Nh (at the source redshift) and of the energy selection band. 
Therefore, a given source with an intrinsic count rate above the 
survey threshold could have been missed in the survey since, 
because of the absorption, its observed count rate falls below 
the survey count rate threshold. 

In order to compute the Zma.K of each source (absorbed or 
not) in section 3 we have used the observed MOS2 count rate in 
the selection band (4.5-7.5 keV) and the count rate limit of the 
HBSS survey in the same energy band (= 2x10"^ cts/s). The in- 
trinsic (2-10 keV) luminosities (as derived from the X-ray spec- 
tral analysis of each source) have been thus used to divide the 
sources into luminosity bins and to derive the luminosity func- 
tion. A similar procedure was recently used bv Shinozaki et akl 
in computing the local luminosity function of the ab- 
sorbed an d unabsorbed AGN in the HEAOl surveys; however 
Shinozaki et al. (2006) did not prove the correctness of this ap- 
proach and, in particular, if this approach may recover the true 
space density of the absorbed AGN population. 

To investigate this aspect we have simulated a source sam- 
ple assuming: 

a) a 2-10 keV local luminosity function described by a 
smoothly connected two power-law function (see section 3.2) 
having the following parameters yi = 1.55, 72 - 2.61, 
LogL* - 44. These parameters are those describing the in- 
trinsic 2-10 keV XLF of the absorbed AGN population in the 
HBSS survey; 

b) a LDDE cosmological evolution model (see section 3.1) 
having best fit parameters as follows: pl=6.5, p2=-1.15, z* = 
2.49, a=0.20. Log L„=45.80; 

c) a flat intrinsic Log Nh distribution from 10^*^ to 10^^ 

-2 

cm . 

Folding all this information together we have produced a 
simulated sample of ~ 40000 sources that would be visible at 
the count rate threshold of the HBSS survey. For each source 
we have computed the observed count rate in the 4.5-7.5 keV, 
the redshift, the absorbing column density Nh and the intrin- 
sic luminosity in the 2-10 keV band. For simplicity we have 
assumed the same intrinsic photon index (F = 1 .9) for all the 
sources. This simulated sample has been thus analyzed with the 
same procedure used for the real HBSS source sample, namely, 
using the observed 4.5-7.5 count rate (and the count rate limit 
of the HBSS survey), combined with the input LDDE model, 
to evaluate space densities and the intrinsic (2-10 keV) lumi- 
nosities to split the sources into luminosity bins. 




43 44 45 



Log (2-10 keV; erg s-') 



Q 

cu 
o 




0) 



i. 0.6 - 



43 44 45 

Log L^ (2-10 keV; erg s-') 

Fig. A.l The results reported in this Figure are those relative 
to the sources filtered by an absorbing column density Nh be- 
tween lO^-' and 10^"* cm"^. Upper panel: comparison between 
the input XLF (solid line) and the XLF computed (dots) using 
the method discussed here. Lower panel: ratio between the in- 
put and output XLF. The sohd lines enclose Icr error bars on 
this ratio. 

The results obtained for one extreme example (the sources 
having an absorbing column density in the range between 10^^ 
and 10^'* cm"^) are shown in Figure A.l; we have reported the 
input XLF, the output XLF and the ratio between the derived 
XLF at z=0 from the simulated sample and the input local XLF. 
Is is worth noting the very good agreement between the input 
and output XLF. Similar results have been obtained for less ex- 
treme values of the Nh- For the purpose of this paper the main 
conclusion we reached using these simulations is that for the 
HBSS selection band (4.5-7.5 keV) our approach to compute 
space densities can recover the true space density of the AGN 
population as a function of the intrinsic Nh up to Nh ~ 10^^ 
cm"^, even if in the A^^ bin between lO^^ and lO^"* about 50 % 
of the AGN are lost (i.e. observed count rate is below the survey 
threshold) because of the intrinsic absorption. Our simulations 
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also show that the above statement is not true for A^^ above 10 
where a dramatic number of objects (~ 95% between 10^* and 
10^^) are lost because of absorption (note that for these high A^^ 
values the situation could be even worst if Compton scattering 
is considered). In summary using the HBSS survey and our 
computation method we can recover the true space density of 
absorbed AGN up to Nh ~ cm"^, while we are quite inef- 
ficient in selecting AGN with an intrinsic Nh above this value. 
The shape of the intrinsic XLF, at least for AGN with Nh up 
to ~ 10^'* cm"-^, is therefore recovered allowing us to compare 
the XLF of absorbed and unabsorbed AGN consistently. We 
have repeated these simulations using different values of the 
parameters for the XLF and/or for the cosmological evolution 
obtaining the same results thus giving us confidence about the 
overall procedure. 



